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ABSTRACT 

A mechanism based on the penetration of interstellar ultraviolet photons into the inner layers of 
clumpy circumstellar envelopes around AGB stars is proposed to explain the non-equilibrium chem- 
istry observed in such objects. We show through a simple modelling approach that in circumstellar 
envelopes with a certain degree of dumpiness or with moderately low mass loss rates (a few 1CU 7 
M© yr _1 ) a photochemistry can take place in the warm and dense inner layers inducing important 
changes in the chemical composition. In carbon-rich objects water vapor and ammonia would be 
formed with abundances of 10 -8 — 10~ 6 relative to H2, while in oxygen-rich envelopes ammonia and 
carbon-bearing molecules such as HCN and CS would form with abundances of 10 -9 — 10~ 7 relative 
to H2. The proposed mechanism would explain the recent observation of warm water vapor in the 
carbon-rich envelope IRC +10216 with the Herschel Space Observatory, and predict that H2O should 
be detectable in other carbon-rich objects. 

Subject headings: astrochemistry — molecular processes — stars: AGB and post- AGB — circumstellar 
matter 



1. INTRODUCTION 

Evolved stars in the asymptotic giant branch (AGB) 
phase undergo important mass loss processes that pro- 
duce extended circumstellar envelopes (CSEs) composed 
of dust and molecules. The molecular composition of 
CSEs is affected by several processes during the travel 
towards the interstellar medium (ISM) but it is originally 
established in the stellar atmosphere under thermochem- 
ical equilibrium (TE) conditions and is to a large extent 
governed by the C/O abundance ratio. In C-rich stars 
(C/O > 1) most of the oxygen is in the form of CO re- 
sulting in a carbon-based chemistry while in O-rich stars 
(C/O < 1) CO l ocks most of the carbon leaving little for 
other molecules (|Tsuiilll973t see also Fig. [l}. 

Astronomical observations have in the main confirmed 
this picture, although with a significant number of dis- 
crepancies. HCN emission is widely observed in O- 
rich objects dBujarrabal et al.lH99l iBieging et~aT1 120001: 
iDecin et al. 2008), i n some of them coming fr om the inner 
regions of the CSE (jDuari fc Hatchellll2000h . Water va- 
por in the C-rich envelope IRC +10216 is confined to the 
warm near surroundin gs of the star, acc ording to recent 
Herschel observations ( Decin et al.ll2010T ). Moreover, am- 
monia is observed in the inner regions of C- and O-rich 
CSEs (jKeadv fe Ridgwavlll993llMenten fe Alcolealll995H 
with abundances much larger than predicted by TE. 

On the theoretical side, shocks induced by the stel- 
lar pulsation have been proposed to explain the non- 
equilib rium chemistry o bserved in the inner regions of 
CSEs (Chcrchncfi 2006). These models explain the for- 
mation of HCN and CS in the inner wind of O-rich 
CSEs, although fail to explain the presence of H2O in 
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IRC +10216 and of NH 3 in both C- and O-rich CSEs. 

In this Letter we investigate an alternative mechanism 
of non-equilibrium chemistry, based on the penetration 
of interstellar ultraviolet (UV) photons into the inner 
regions of CSEs with a certain degree of dumpiness. 

2. MODEL 

The model is based on a central AGB star surrounded 
by a spherical circumstellar envelope, whose physical pa- 
rameters are given in Table [TJ and is adopted to inves- 
tigate the chemistry in both C- and O-rich CSEs with 
mass loss rates between 10 -7 and 10 -5 M© yr _1 . 

In CSEs with an intense mass loss process strictly 
isotropic and homogeneous, the inner regions are 
well shielded from interstellar UV light, so that only 
the outer layers are affected by photochem istry (e.g. 
ICharnlev et all fl995l: IWillacv fc Mihan [1997m . Obser- 
vations, however, have shown that CSEs have usu- 
ally clumpy structures both at small and large scales 
T99 l IGuelin et al.llT997l: IWeigelt et al.l 
200a ILeao et al.ll2006fl . which allow for 



(Chapman et al. 
1998: Fong et al 



a deeper penetration of interstellar UV photons into the 
inner regions. To model the effects of dumpiness on the 
circumstellar chemistry we adopt a simple approach in 
which the CSE consists of two components: a major one 
whose inner regions are well shielded against interstellar 
UV light, and a minor one (which accounts for a frac- 
tion f m of the total circumstellar mass) for which the 
shielding matter located in the radial outward direction 
is grouped into clumpy structures leaving a fraction 
of the solid angle of arrival of interstellar photons free of 
matter. 

The gas phase chemical composition of these two com- 
ponents is computed as they expand from the inner- 
most regions (r = 2i?„) until the end of the CSE. The 
adopted abundances at the initial radius are given in 
Table [5] 13 CO is also i ncluded with an a bundance 30 
times lower than 12 CO (jMilam et al.ll2009D . The chem- 
ical network has been used in previous chemical models 
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TABLE 1 

Model physical parameters 



Star effective temperature (T») 
Stellar radius (R*) 
Expansion velocity (v) a 

Gas kinetic temperature (T) b 
Volume density of gas particles (n) c 



5 km s 
15 km s 



2000 K 
5 X 10 13 cm 
1 for r < 5 R* 
- 1 for r > 5 R* 
n(r/R t )-°- 7 
M / (A-nr 2 (m s )v) 



Note. — r is the radius measured from the center of the star. 
a Velocity field similar to that adopted in previous studies of in- 
ner CSEs dKeadv fc RidgwavlH993l : IFonfria et al.H20Qgj). b Values of 
the exponent are typically between —0.5 and —1 (Justtanont et al. 
ll994tlFon"fn'a et al.ll200l) . " M is the mass loss rate and (m g ) is the 
mean mass of the gas particles (typically ~ 2.3 times the mass of the 
hydrogen atom in CSEs). 



of warm and dense UV illu minated regions ([Cernicharol 
l2004IAgundez et al]|2008al) . Photodissociation and pho- 
toionization rat es are evaluated as a function of the visual 
extinction Ay (jWoodallet all [20071: Ivan Dishoeck et al.l 
I2006D , adopting the interstellar UV field of lDrainei (|1978D . 
For 12 CO the photod i ssocia tion rate is evaluated accord- 
ing to IMamon et al.l ((1988) , who included the effect of 
self-shielding, and for 13 CO through the expr ession 2 x 
10- 10 exp(-2.5^ y ) s" 1 (|Woodall et al.ll2007T) . 

For the major component, shielded by a smooth enve- 
lope, Ay depends on the column density of H nuclei Ah 
in the r adial outward direc tion as Ay — Nn(cm^ 2 )/1.87 
x 10 21 (|Bohlin et alj|l978l ). The UV field for this com- 
ponent may be expressed as: 



(i) 



where I\ is the unattenuated interstellar UV field at a 
wavelength A, fl the solid angle of arrival of most of the 
UV flux (depends strongly on the radius and ranges from 
a small fraction of ir sr in the inner regions up to almost 
47r sr in the outermost layers), and [yU/Ay] the ratio of 
the dust extinction at A and at visual wavelengths (3.8 for 
A = 1 250 Aas found for the ISM bv lFitzpatrick fc Massal 
1990). For the minor component, shielded by a clumpy 
envelope, the UV field may be analogously expressed as: 



AttJ x = i x n 



(l-.fc)cxp{-^^i^}+/o] (2) 



For the minor component we compute an effective visual 
extinction Ay by substituting Ay by Ay into Eq. (TTJ 
and equalling the right parts of Eqs. ([l} and ((2|). 

Near-infrared interferometry, able to trace the dust 
emission at milli-arcsecond scales, has revealed a ex- 
tremely clumpy structure in objects such as IRC +10216, 
with five individual clumps in the close surroundings of 
the star, some of them with angle separations of up to 
20-30° (jWeigelt et al.lll998| ). Millimeter- wave interfer- 
ometry of molecular lines can, unlike infrared observa- 
tions of dust, provide information on the projected ve- 
locity in the line of sight and allow to build three di- 
mensional maps. Observations of IRC +10 216 in molec- 
ular lines of CN, C?H, C 4 H, an d HC 5 N (|Guelin et al.l 
119931: iDinh-V-Trung fc Limll2008l ) have shown that these 
species are distributed in a circumstellar shell with a ra- 
dius of 15-20", and with two conical holes in the NNE 
and SSW directions having an aperture angle of about 
30°. This corresponds to a solid angle of tt/A sr, which 



TABLE 2 

Initial abundances relative to H2 in C- and O-rich CSEs 



Carbon-rich 



Oxygen-rich 



Species 


Abundance 


Rei 


Species 


Abundance 


Rei 


He 


0.17 




He 


0.17 




CO 


8 X 10~ 4 


(1) 


CO 


3 X 10~ 4 


(1) 


N 2 


4 X 10~ 5 


(2) 


N 2 


4 X 10~ 5 


(2) 


C2H2 


8 x 10~ 5 


(3) 


H 2 


3 x 10~ 4 


(-1) 


HCN 


2 x 10~ 5 


(3) 


co 2 


3 x 10~ 7 


(5) 


SiO 


1.2 x 10~ 7 


(3) 


SiO 


1.7 x 10~ 7 


(6) 


SiS 


10~ 6 


(3) 


SiS 


2.7 x 10~ 7 


(7) 


CS 


5 x 10~ 7 


(3) 


so 


lO" 6 


(8) 


SiC 2 


5 x 10" 8 


(3) 


H 2 S 


7 x 10" 8 


(9) 


HCP 


2.5 x 10" 8 


(3) 


PO 


9 x 10~ 8 


(10) 


References. — (1) ITevssier et al. 


j2006tl; (2) TE abundance 


; (3) 



JSOOSft; (5) ITsuii et~aTl 1119971): (6) ISchoier et all 112009); ~(7) 
Schoier et al. (2007); (8) Bujarrabal et al. (1994): (9) Ziurys ct al. 
120071) ; (lO) ITenenbaum et al.l 120071) . 

may be a good fraction of fl in the inner circumstellar 
regions. Both f m and fa are phenomenological parame- 
ters, difficult to quantify in any CSE. Anyway, adopting 
fm = 0.1—0.2 (the minor UV illuminated component ac- 
counting for just a 10—20 % of the total circumstellar 
mass) and fa = 0.2—0.5 (in line with the observational 
data described above for IRC +10216) seems reasonable 
for a CSE with a sufficient degree of dumpiness and al- 
lows for a sufficient penetration of interstellar UV pho- 
tons into the inner layers of the CSE. For all the models 
we have adopted f m — 0.1 and fa = 0.25. 

3. RESULTS 

Fig. [1] shows the TE abundance of various molecules 
calculated under conditions typical of the stellar at- 
mosphere of AGB stars. The ca lculations have been 
done with the code described in iTeiero fc Cernicharol 
(|1991|) and most of the thermochemical data have been 
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Fig. 1. — TE abundances as a function of temperature for a 
C- and O-rich gas (with a C/O abundance ratio of 1.5 and 0.5 
respectively) with a constant particle density of 10 14 cm -3 . 



Photochemistry in the inner layers of clumpy circumstellar envelopes 



10 M Q yr 




10 10 
radius (cm) 

Fig. 2. — Calculated abundance of several molecules as a function of radius for carbon-rich CSEs with mass loss rates of 10 — , 10 — 6 , and 
10~ 7 Mq yr -1 . Dashcd-dotted lines correspond to the abundance of the minor UV exposed component properly corrected by the factor 
/m, while continuous lines correspond to the abundance weighted-averaged over the minor and major components, which can be expressed 
as Xi(r) = (1 - fm.)xf^ (r) + f m x™(r), where xf 1 (r) and x™(r) are the abundances of the species i at a radius r in the major and minor 
components, respectively. Although in the reality the minor UV exposed and the major UV shielded components would not coexist in 
space, this treatment allows to simplify taking the abundance averaged over all radial directions. Note that in those regions where the 
dashed-dotted and continuous lines of a given species coincide, only the minor UV exposed component contributes to the abundance. 



taken from |Chase] dl99j$. The main purpose of these 
calculations is to show that several molecules which are 
observed in the inner regions of CSEs have very low TE 
abundances (e.g. H2O, NH3, and SiH4 in C-rich objects, 
and HCN, CS, and NH 3 in O-rich objects). 

Now focusing on the models based on chemical kinetics, 
Fig. [2] and Fig. [3] show the calculated abundance distri- 
bution of some molecules in C- and O-rich CSEs with 
mass loss rates of 10 -5 , 10 -6 , and 10 -7 M yr -1 . 

In C-rich CSEs with mass loss rates as high as 10~ 5 
M Q yr -1 , typical of an object such as IRC +10216, wa- 
ter vapor would be effectively formed in the dense and 
warm inner regions of the minor UV illuminated compo- 
nent, for which Ay is < 1 mag, with a global abundance 
relative to H 2 in excess of 10~ 7 (see Fig. [2]). In such re- 
gions the photodissociation of 13 CO and SiO, the major 
reservoirs of oxygen besides 12 CO (hard to photodisso- 
ciate due to self-shielding effects), liberates atomic oxy- 
gen which is effectively converted into water through the 
chemical reactions 



O + H 2 ->• OH + H 
OH + H 2 -» H 2 + H 



(3) 
(4) 



which despite having activation barriers are rapid enough 
due to the high temperatures attained in these inner 
layers. The photodestruction of molecules in the minor 
UV illuminated component occurs fast, but for some of 
them (e.g. H 2 0) the formation rate is high enough to 
allow them to extend up to relatively large radii, ~ 10 16 
cm (see Figj2]). This mechanism would explain the ob- 
servation with Hcrschel of warm water vapor in the in- 
ner circumstellar r egions of the carbon star IRC +10216 
(jDecin et al.l I2010T ). Other mech anisms proposed suc h 
as sublimation of cometary ices (|Melnick et al.l[200l . 
Fischer- Tropsc h catalysis on the surface on iron grains 
or radiative asso ciation between O and 
Cernicharo 2006) , place water in cool re- 



gions located farther than 10 15 cm. 

The proposed mechanism would also have some other 
interesting consequences. Hydrides other than H 2 0, 
such as NH 3 , CH 4 , H 2 S, SiH 4 , and PH 3 , could also 
be effectively formed in the inner CSE by successive 
hydrogenation reactions of th e heavy atom. All them 
are observed in IRC +10216 dKeadv fc Ridgwavi Il99l 
ICernicharo et al.ll2000t lAgundez et al Il2008b[ ) with abun- 
dances which, except for CH4, are much larger than pre- 
dicted by TE (see Fig. [1]). Ammonia, for example, is 
observed in the inner CSE of IRC +1 0216 with an abun- 
danc e relative to H 2 of 10~ 7 — 10~ 6 (|Keadv &; Ridgwavi 
Il993t lHasegawa et al.l I2006D . and yet no efficient for- 
mation mechanism has been proposed, apart from the 
suggestion that it could be formed on grain surfaces 
(jKeadv fe Ridgwavi 1 19931 ). Our model predicts an effec- 
tive formation for NH 3 and H 2 S (see Fig. [2]), but not for 
SiH4 and PH3 (likely due to the lack of chemical kinet- 
ics data for the relevant hydrogenation reactions) . Other 
molecules such as HC3N increase their abundance in the 
inner envelope due to the penetration of interstellar UV 
photons (see Fig. [2]), something that has been recently 
confirmed through observations of IRC +10216 at A = 0.9 
mm with the IRAM 30-m telescope (jDecin et al.l 120101: 
Kahane et al. in preparation). 

Still focusing on C-rich sources (see Fig. [2]), if we move 
toward lower mass loss rates then the whole CSE starts 
to be more transparent to interstellar UV photons. For 
example, for a mass loss rate as low as 10 -7 M Q yr -1 
most of H 2 is formed in the major UV shielded com- 
ponent, which is no longer shielded as it has a visual 
extinction < 1 mag. Thus, for C-rich CSEs with mod- 
erately low mass loss rates (up to a few 10 -7 M yr -1 ) 
we should expect a relatively large H 2 abundance even 
if the CSE is not particularly clumpy, prediction that 
should be easily tested with Herschel. 

In the case of O-rich CSEs the penetration of interstel- 
lar UV photons into the inner layers has also interesting 
chemical effects (see Fig. [3]). Among them it is worth 



Agiindez et al. 



10 7 M yr 1 




10 10 



10 10 
radius (cm) 



Fig. 3. — Same as in Fig.[2]but for oxygen-rich CSEs. 



mentioning the formation of NH3, CH4, HCN, and CS 
in the inner envelope with abundances relative to H2 in 
the range 10 -8 — 10~ 7 , i.e. much larger than predicted by 
TE calculations (see Fig.[I|. The formation of C-bearing 
molecules in a dense a nd warm UV illuminate d O-rich gas 
has been discussed bv lAgundez et al.l (|2008al ) in the con- 
text of the chemistry of protoplanetary disks. HCN, CS, 
and NH3 are observed in O-rich CSEs with abundances 
of lQ- 7 -lQ- 6 relative to H? (iBuiarrabal et aLI 119941: 
iMenten fc Alcoleal fl995t iBieging et all 120001 ). which are 
somewhat higher than predicted by us. Other mecha- 
nisms based on sh o cks induced by t he stellar pulsation 
(jDuari et al.l 119991: iCherchneffl E006) predict fractional 
abundances for HCN and CS of 10~ 6 — 10~ 5 , which are 
somewhat higher than observed, but a negligible abun- 
dance for NH 3 . 

4. CONCLUSIONS 

We have shown through a simple modelling approach 
that in CSEs envelopes with a certain degree of dumpi- 



ness or with moderately low mass loss rates (a few 10 7 
Mq yr -1 ) a photochemistry can take place in the warm 
and dense inner layers inducing important changes in 
the chemical composition. This mechanism allows for 
the formation of H2O and NH3 in C-rich objects and 
HCN, CS, and NH 3 in O-rich objects, with abundances 
much higher than predicted by thermochemical equilib- 
rium but close to the values typically derived from astro- 
nomical observations. This mechanism explains the re- 
cent observation of warm water vapor in the carbon-rich 
envelope IRC +10216 with the Herschel Space Observa- 
tory, and predict that H2O should be detectable in other 
carbon-rich objects. 
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